Metal-insulator quantum critical point beneath the high T c superconducting dome 
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An enduring question in correlated systems concerns whether superconductivity is favoured at 
a quantum critical point (QCP) characterised by a divergent quasiparticle effective mass. Despite 
such a scenario being widely postulated in high T c cuprates and invoked to explain non-Fermi 
liquid transport signatures, experimental evidence is lacking for a critical divergence under the 
superconducting dome. We use ultra- strong magnetic fields to measure quantum oscillations in 
underdoped YBa2Cus06+x, revealing a dramatic doping-dependent upturn in quasiparticle effective 
mass at a critical metal-insulator transition beneath the superconducting dome. Given the location 
of this QCP under a plateau in T c in addition to a postulated QCP at optimal doping, we discuss 
the intriguing possibility of two intersecting superconducting subdomes, each centred at a critical 
Fermi surface instability. 

PACS numbers: 71.45.Lr, 71.20.Ps, 71.18.+y 



A continuous zero temperature instability between 
different groundstates— termed as a quantum critical 
point— is characterised by a divergence in a relevant sus- 
ceptibility [1]. In strongly correlated systems @], the in- 
fluence of criticality on the entire body of itinerant elec- 
trons results in a global divergence of the effective mass— 
which is recognized as the key defining experimental sig- 
nature of quantum criticality @, S] • The growth of elec- 
tronic correlations on the zero temperature approach to 
the critical instability can be experimentally accessed by 
the tuning of parameters such as pressure and doping. 
Quantum oscillation measurements are ideally suited to 
investigate the effects of such tuning due to the direct 
access they provide to the effective mass of the elemen- 
tary fermionic excitations that can be traced across the 
Quantum Critical Point (QCP) Such a direct probe 
is crucial in superconducting materials, where bulk ther- 
modynamic signatures of quantum critical behaviour of 
the normal quasiparticles [l|, [2|, [5[ are difficult to access 
due to the overlying superconducting dome. 

While the emergence of high T c superconductivity in 
the cuprate family is inextricably linked to the parent 
Mott insulating compound, remarkably little is known 
about the physics of the metal-insulator crossover 
and its relation to electronic correlations. By using 
quantum oscillation measurements in strong magnetic 
fields to access normal state quasiparticles in under- 
doped YBa2Cu306+^, we uncover a striking doping- 
dependent upturn in the effective mass at the loca- 
tion of the metal-insulator crossover 0, [H, 0, EE 11]. 
Our findings provide bulk thermodynamic evidence for 
a metal-insulator quantum critical poin t (QCP) in high 
T c cuprates [H, EE EI EE EE EE Ell, without requir- 
ing extrapolation below the superconducting dome. The 
effective mass divergence unaccompanied by a change in 



Fermi surface area away from half-filling signals a novel 
many-body mechanism [2(j that drives insulating be- 
haviour in underdoped cuprates. 

We trace the doping dependence of quantum oscil- 
lations with increased underdoping of YBa2Cu306+cc 
(x=0.54, 0.51, 0.50, 0.49). Of the multiple Fermi sur- 
face orbits detected in a subset of samples of x =0.50, 

(see Ap- 



0.51, and 0.54 compositions [2l|, [22|, [23|, [24|, |2£ 
pendix Fig. [8] for an example of the higher /3-frequency 
observed in current measurements on the x =0.54 dop- 
ing), we focus on the a pocket of carriers that shows 
the most prominent quantum oscillations in all measured 
compositions. Figure [T] shows examples of quantum os- 
cillations we measure using contactless methods, where 
changes in the resistivity are reflected as a shift in reso- 
nance frequency (A/) of an oscillator circuit. Measure- 
ments are made down to temperatures of 1 K using a 
tunnel-diode oscillator in a slowly swept magnet reaching 
fields of 55.5 T, and a proximity detector in a two-stage 
magnet system reaching fields of 85 T (see Appendix). 
The high magnetic fields used here enable access to the 
evolution of low energy quasiparticle excitations by the 
suppression of superconductivity (wherein we refer to the 
zero resistance state); the crossover field into the high 
magnetic field resistive state (H r ) is shown as a function 
of oxygen composition (x) and temperature (Fig. 1 in- 
set). The location of H r is close to the irreversibility field 
H 1VV determined by torque measurements in our previous 
work (22], 

The quasiparticle effective masses (m*) are extracted 
by performing a Lifshitz-Kosevich fit to the temperature 
dependent amplitude of the observed quantum oscilla- 
tions [26[ (shown in Fig. [2]). Our key experimental find- 
ing is that m* exhibits a steep upturn in samples of pro- 
gressively lower oxygen concentration x (Fig. [3]). The 
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FIG. 1: An example of the measured magnetic field- 
dependent resonance frequency change Af (see text) for one 
sample of each composition at 1.5 — 1.6 K measured in both 
55.5 T (shown by solid lines) and 85 T (shown by dotted 
lines) magnets (see Appendix). The corresponding sample 
(i.e., A, D, G, H, I) is indicated in each case. Single crys- 
talline platelets of average dimensions ~ 800 x 500 x 50 fim 
are coupled inductively to the face of a coil of ^ 5 turns (with 
its axis parallel to H) that forms part of the contactless con- 
ductivity circuit (see Appendix). H causes a crossover into a 
high magnetic field resistive state whose in-plane skin depth 
(rsj 100 /im at 46 MHz) increases the coil inductance causing / 
to drop. The resistive crossover field H T (also plotted versus T 
using solid symbols in the inset) is determined from the max- 
imum in the derivative d 2 f/dB 2 . Hollow circles show static 
field measurements of H r measured on a sample of one of the 
same compositions (x = 0.54) in a 45 T static field magnet. 
Given the thermal energy generated by flux flow motion in 
Type II superconductors, we ensure sample immersion in the 
liquid cryogen to ensure heat dissipation and well controlled 
temperatures (see Appendix). 

masses are independent within fit uncertainties of sam- 
ple (different crystals of the same composition), magnet 
system (sweep rate), magnetic field range, distance from 
the irreversibility field, and experimental setup (see Ap- 
pendix). The decrease in quantum oscillation amplitude 
with deoxygenation beyond that expected for the increase 
in mass indicates a Landau level broadening— associ- 
ated either with increased oxygen disorder, a stronger 
pairing potential or an increased probability of scatter- 
ing reflecting the increase in m*. In contrast to the 
striking increase in m*, the cross-sectional area of the 
pocket Ak = (27re /fi)F (where F is the observed quan- 
tum oscillation frequency in reciprocal magnetic field 
1/B « 1/fjLoH 0) shows a comparatively weak depen- 
dence on x (see Appendix). 

A tuning-driven divergence in m* is identified by a 
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FIG. 2: Quantum oscillations are measured on nine de- 
twinned ortho-II ordered YBa2Cus06+x single crystals using 
the contactless conductivity technique in the motor generator- 
driven 55.5 T magnet and 85 T magnet, with nominal com- 
positions x = 0.49 (samples A, B and C), 0.50 (samples D, 
and E), 0.51 (samples F and G) and 0.54 (samples H, I and 
J). The figure shows fits of the Lifshitz-Kosevich expression 
(a = aoX/ sinhX , where X = 27T 2 m*kBT/heB [26|) made 
to FFT amplitudes of quantum oscillations in seven samples 
(i.e., A, B, D, E, F, G, H and I) measured in the 55.5 T (a) 
and 85 T (b) magnet over the field ranges shown in Fig. [3] 
as a function of temperature (T). Amplitudes for each of the 
dopings are renormalized to coincide at 4 K for comparison 
purposes. Error bars correspond to the noise floor of the FFT. 
The inset shows examples of the measured oscillations at sev- 
eral different temperatures for sample 6.49 (A) measured in 
the 55.5 T magnet, and samples 6.49(A) and 6.50(D) in the 
85 T magnet. The same analysis performed over a subset of 
field ranges is shown in the Appendix. 



collapse in the inverse many-body mass enhancement 
(mb/m*, where is the band mass), and hence Fermi 
temperature (Tp) to zero at a critical value of the tun- 
ing parameter [3] . Figure 3d shows the ratio / m* as a 
function of x in YBa2Cu306+^ (rrih ~ 0.5 m e is estimated 
from conventional band theory [22] and is assumed to re- 
main constant for the incremental changes in x accessed, 
given the largely unchanged pocket area), and Fig. [H 
shows the inferred Fermi temperature Tp = TieF jm*k^. 
A precipitous linear drop in these quantities is seen with 
reduced oxygen concentration, presaging their vanishing 
in the vicinity of a critical doping x c . Linear interpola- 
tion yields x c « 0.46 as the location of a putative quan- 
tum critical point beneath the superconducting dome in 
YBa2Cu306+x (seen from Fig. [3]i). 

Critical behaviour tuned by doping rather than by 
magnetic field is evidenced by the absence of a discernible 
magnetic field dependence of either F or m* in the range 
fioH ~26 to 85 T. Thermal conductivity measurements 
in zero field on YBa2Cu306+^, also show a notable drop 
for x < x c [7| (see Appendix), and /isr measurements 
at zero field reveal an abrupt change in /isr line shape 
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FIG. 3: Examples of quantum oscillations in each of the sam- 
ples in Fig. [T]at T ^ 1.5 — 1.6 K measured in the 55.5 T (a) 
and 85 T (b) magnets after background polynomial subtrac- 
tion. The largest amplitude quantum oscillations of frequency 
F a ~ 535 T are observed consistently in all compositions, 
c, Effective mass m* of the quantum oscillations (where m e 
refers to the free electron mass) extracted from fits shown 
in Fig. [2] (blue and pink circles referring to the 55.5 and 
85 T magnets respectively) plotted as a function of x. d, 
Inverse many-body mass enhancement (rrih/m*) and a linear 
fit shown on the left-hand axis as a function of x. The super- 
conducting temperature T c as a function of x shown by dots 
(see Appendix) together with measured T c values for sam- 
ples A through J (using a SQUID magnetometer) shown by 
diamonds on the right-hand axis. 
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FIG. 4: Diagram in which the Fermi temperature Tf is ex- 
tracted from the measured values of m* and F (see text), with 
blue and pink circles referring to data taken in the 55.5 T 
and 85 T magnets respectively. Green diamonds represent 



the metal-insulator transition temperature Tmi as described 
in the text (extracted from Refs. |a,E|)- Fits to Tf and Tmi 
are represented by dotted lines and a sharp demarcation be- 
tween different coloured regions. The intercepts x c and xmi 
refer to the extrapolations of Tf and Tmi to T = 0, indicating 
the existence of a putative quantum critical point. The solid 
red line depicts twice the superconducting temperature T c of 
YBa2Cu306+cc at zero field (from the Appendix), while the 
dark grey line represents its extrapolation to notional higher 
dopings (based on an analogy with S^-ccLaxCuCU)- The 
fuzzy conical region centred at optimal doping represents an- 
other quantum critical region postulated to occur where the 
Fermi surface consisting of small pockets transforms to a large 
Fermi surface (see text). The red dashed curves represent a 
possible division of the superconducting dome into two inter- 
secting subdomes overlying two distinct critical points: x c , 
and a postulated critical point at optimal doping. 



below x c Intriguingly, the critical doping x c is lo- 
cated at the same region of doping where the postulated 
crossover from metallic to insulating behaviour of the 
normal carriers 0, S, 0, OH EH onsets, characterized by 
a low temperature logarithmic divergence in the resis- 
tivity [1, EH EH • We extract in Figure [H the doping- 
dependence (for x < x c ) of the metal-insulator transi- 
tion (or crossover) temperature (TMi)from the in-plane 
resistivity data reported in Refs. [3, Here Tmi de- 
notes the temperature at which the zero (high) field re- 
sistivity reaches its lowest value before logarithmically 
diverging at low temperatures. From Fig. 2J we find 
that with increasing doping, Tmi collapses linearly to- 
ward xmi ~ 0.46, which denotes the zero temperature 
metal-insulator transition. The coincidence of xmi an d 
x c at T = signals an association of the experimentally 
observed collapse in Fermi temperature and divergent ef- 
fective mass with a zero (low) temperature continuous 
metal-insulator transition. The lack of saturation in the 
logarithmically diverging resistivity for x < x c 0, EH 11] 
indicates that the metal-insulator QCP demarcates a 
sharp transformation of the entire body of conduction 
electrons from small Fermi surface pocket to insulating 



regime in YBa2Cu306+^- 

Although it was previously considered that disor- 
der (e.g. weak localisation) or band depopulation 
drives insulating behaviour at x < x c in underdoped 
YBa2Cu306+£ [l9[, the steep upturn we observe in m* 
unaccompanied by a change in F signals that elec- 
tron correlations are central in driving the develop- 
ment of insulating behaviour for x < x c in underdoped 
YBa2Cu306+^ [20[. While a correlation-driven metal- 
insulator transition is not entirely surprising given the 
proximity to the Mott insulating regime dominated by 
Coulomb repulsion, the continuous nature of the metal- 
insulator transition (indicated by the collapse of the resis- 
tivity upturn temperature at xmi) is unexpected [6]. Fur- 
thermore, the location of the observed correlation-driven 
metal insulator transition away from half filling (x = 0) 
in YBa2Cu306+;r suggests an alternate theoretical sce- 
nario (e.g. Refs. [27|, [28[) to that originally proposed by 
Brinkman and Rice [20j]. One possibility is the interplay 
of additional interactions other than those considered in 
the Brinkman-Rice picture. Signatures of magnetic or- 
der have been reported on both sides of x c ; a collapse 
in spin excitation gap has been reported as x is reduced 
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below x c [Io[, while spin density wave ordering has also 
been suggested to be responsible for Fermi surface recon- 
struction at x > x c and /iqH > 30 T [29]. A contender 
for an order parameter that onsets below x c to drive the 
continuous metal-insulator transition at finite doping is 
charge order - static charge order has been observed to 
develop below x c by inelastic neutron scattering measure- 
ments [30] . A transformation between local and itinerant 
magnetism near x c may be indicated, the development of 
local magnetic moments below x c having been reported 
from jisr experiments [8|- 

While low temperatures are required for our measure- 
ment of quantum oscillations— with their observation re- 
quiring k^T to fall well within the Landau level spac- 
ing (h(j c /k B ~ 80 K at B = 85 T and x = 0.54)- 
the Fermi energy scale (Tp) associated with the observed 
Fermi surface pockets extends to energies greatly exceed- 
ing T c away from the QCP. The rapid collapse of this 
energy scale at x c — located under the local maximum 
(or plateau) of the superconducting dome— mirrors the 
behaviour seen in strongly correlated f-electron supercon- 
ductors, in which case a diverging effective mass has been 
reported at a QCP |5[ under the superconducting dome 
maximum, where f-electrons are removed from participa- 
tion in the Fermi surface volume [4(. Finally, we note 
that another quantum critical point (or extended region 
of criticalit y) i s postulated to occur near optimal dop- 
ing [H 0,13, EE El, E3, E3 (i.e. the maximum of the 
upper superconducting dome in YE^CusOe+cc) where 
the small Fermi surface pockets in YE^CuaOe+a, [2lll22| 
are expected to evolve into a large Fermi surface recently 
observed in Tl2Ba2Cu06+5 [31]. An intriguing possi- 
bility to consider therefore is the existence of two in- 
tersecting superconducting domes in high T c cuprates— 
perhaps similar to the seminal heavy fermion supercon- 
ductor CeCu2Si2 [2|— where each of the superconducting 
subdomes is centred at a distinct critical Fermi surface 
instability. 

APPENDIX 

T c and sample compositions: The T c curves plotted in 
Figs. 3d and 4 are taken from Ref. [32[. Non oxygen- 
ordered samples such as those measured in Ref. [7] were 
previously reported to have a slightly different T c versus 
x' dependence to those in Ref. [32[ , with a putative metal- 
insulator transition reported to occur at x' Ml ~ 0.55. In 
this work, for accurate comparison T c values are used as 
a means of renormalizing doping values (x) of ortho-II 
ordered and oxygen disordered samples grown by differ- 
ent methods (following Li et al. in Ref. [lQ[). Using 
renormalized dopings to compare Ref. [7], then x' Ml is 
equivalent to xmi ~ 0.47 for the current samples— close 
to the extrapolated value shown in Fig. 4. 

Magnet systems used in the experiments: Two differ- 
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FIG. 5: Pulse waveform for each of the 55.5 T and 85 T 
magnets used. The slow sweep rate of the 55.5 T magnet 
and up to 36 T in the 85T magnet is due to magnetic field 
generation by a motor generator. The remaining field in the 
85T magnet is provided by a capacitor bank. 

ent magnet systems are used to perform the experiments. 
Experiments extending to 55.5 T in magnetic field are 
conducted in a motor-generator driven magnet with a 
slower sweep rate and longer pulse length (magnetic field 
versus time profile shown in Fig.[5K) than capacitor bank- 
driven pulsed magnets. Contactless conductivity mea- 
surements performed in this magnet use a tunnel diode 
oscillator circuit with a resonance frequency of ~ 46 
MHz [33]. For the experiments conducted in magnetic 
fields extending to 85 T, an"outsert" magnet powered by 
the motor-generator is swept slowly to ~ 36 T, with the 
remaining magnetic field provided by a capacitor bank- 
driven "insert" magnet (magnetic field versus time profile 
shown in Fig.03). For the experiments performed in this 
magnet, the contactless conductivity measurements use 
a proximity detector circuit resonating at ~ 22 MHz [341 ]. 

A slow ramp rate of the magnetic field is important to 
reduce the effects of flux dissipation heating. While ex- 
periments up to 55.5 T retain a slow ramp rate through- 
out the pulse, flux dissipation heating in experiments up 
to 85 T were minimised due to the slow ramp rate up 
to 37 T in which region the critical current for vortex 
pinning is expected to be largest. For both these mag- 
net systems, different cryostats, measurement probes and 
thermometers were used in addition to different contact- 
less conductivity circuits for detection. 

Temperature control and measurement: The increase in 
H r with decreasing T in the inset to Fig. 1 provides an in- 
situ secondary confirmation of the sample temperature: 
the close correspondence between values for H r extracted 
from static magnetic field measurements and those ex- 
tracted using motor-generator-controlled magnetic fields 
demonstrates that the sample is well-coupled to the liq- 
uid cryogen, minimising heating due to vortex motion 
during these measurements. 

The effects of flux dissipation heating are minimised by 
ensuring immersion of the samples in liquid 4 He through- 
out. For temperatures (measured using a calibrated ther- 
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FIG. 6: FFT amplitudes of quantum oscillations over the lim- 
ited field range 44 - 55.5 T (measured in the 55.5 T magnet) 
and the highest field range 60 - 85 T (measured in the 85 
T magnet) for all dopings as a function of temperature (T). 
Amplitudes of each of the dopings have been renormalized 
as in Fig. 2 to coincide at high temperatures. Error bars 
correspond to the noise floor of the FFT. 



mo meter close to the sample) above 2.17 K, the 4 He is 
re-pressured by back filling with 4 He gas after pumping 
to ensure continued immersion of the sample in the liq- 
uid cryogen during the application of the magnetic field. 
The resistive crossover is seen to be reproduced between 
rising and falling field. 

Extended temperature- dependent amplitudes analysis: 
Temperature dependent amplitudes are shown in Fig. [6] 
over an identical restricted field range 44—55.5 T (mea- 
sured in the 55.5 T magnet) and over the highest field 
range 60—85 T (measured in the 85 T magnet) for all 
dopings x. Similarity with the temperature dependence 
of amplitude extracted over an extended field range in 
Fig. 2 indicates that the measured effective masses shown 
in Figs. 2 and 3, and used to infer the value of Tp in Fig. 
4 are independent (to within the quoted error bar) of the 
magnet system, distance from irreversibility field, and 
magnetic field interval over which it is extracted. 

Doping- dependent frequency analysis: Frequencies cor- 
responding to the a pocket were determined by Fourier 
analysis and by fits of the quantum oscillations in Fig. 
2a,b to A = A cos(27rF a /B + 0)exp(-7/£) . The fre- 
quency can be seen to be largely independent of doping 
in Fig. in contrast to the sharp upturn in effective 
mass seen in Fig. [Tfr. 

The higher f3 quantum oscillation frequency: The 
higher /3 frequency has been observed in a subset of 
measured samples using magnetic torque, contact- 
less conductivity using the Tunnel Diode Oscillator 
(TDO), contactless conductivity using the Proximity 
Detector Oscillator (PDO), and specific heat measure- 
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FIG. 7: (a) Frequency of quantum oscillations corresponding 
to the a pocket (F a ) measured as a function of doping in 
YBa2Cu306+a> (b) Effective mass of the quantum oscillations 
(where m e refers to the free electron mass) extracted from fits 
in Fig. 2 (blue and pink circles referring to the 55.5 and 85 
T magnets respectively) plotted versus doping. 




ments [22, 24, 25] . An example Fourier transform of the 
oscillations showing the (3 frequency Fp ~ 1690 ± 20 T 
from sample I of doping x =0.54 measured in the 85 T 
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FIG. 8: FFT of quantum oscillations measured on a single 
crystal of YBa2Cu3 06.54 using the PDO technique in the 85T 
magnet, revealing the frequency Fp « 1690 ± 20 T, which 
is 3.14(5) times as large and has an amplitude 0.013 times as 
small as the most prominent frequency F a ~ 538 ± 5 T. 



magnet is shown in Fig. [HJ 
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